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Abstract: 

We continue the study of the construction of analytical coefficients of the e-expansion of 
hypergeometric functions and their connection with Feynman diagrams. In this paper, we 
show the following results: 

Theorem A: 

The multiple (inverse) binomial sums 

^ 1 Z 3 

TiTfc^itf -!)••■ s ap (j - 1) , 

3=1 (/) J 

where k = ±1, 5 (j) is a harmonic series, £ a (i) = X^i=i tf> anc ^ c * s an y integer number 
are expressible in terms of Remiddi-Vermaseren functions; 

Theorem B: 

The hypergeometric functions 

pF p -\{A+ae',B+b£,\+Bp-\)z) , p F p -i(A+as, \+A p ; B+be; z) , 

are expressible in terms of the harmonic polylogarithms of Remiddi and Vermaseren with 
coefficients that are ratios of polynomials. 
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1. Introduction 

Feynman diagrams are a primary tool for calculating radiative corrections to any pro- 
cesses within the Standard Model or its extensions. With increasing of accuracy of mea- 
surements, more and more complicated diagrams (with increasing number of loops and 
legs and increasing numbers of variables associated with different particle masses) must be 
evaluated. The essential progress in such calculations is often associated with the inven- 
tion of new (mainly mathematical) algorithms (e.g. Refs. [^, ||]) and their realization as a 
computer programs (e.g. Refs. j|, ||). One fruitful approach to the calculation of Feynman 
diagrams is based on their representation in terms of hypergeometric functions || or mul- 
tiple series |7|, [|. We will refer to such representations as hypergeometric representations 
for Feynman diagrams. Unfortunately, there does not exist a universal hypergeometric 
representation for all types of diagrams. Constructing these representations is still a mat- 



ter of the personal experience of the researcher |9], 10, 11, 12, |13|] . Nevertheless, existing 
experience with Feynman diagrams leads us to expect that all Feynman diagrams should 
be associated with hypergeometric functions. 

For practical applications, finding a hypergeometric representation is not enough. It 
is necessary to construct the so-called e-expansion, which we may understand as the con- 
struction of the analytical coefficients of the Laurent expansion of hypergeometric func- 
tions around rational values of their parameters. In this direction, very limited results are 
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available. 1 The pioneering systematic activity in studying the Laurent series expansion 
of hyper geometric functions at particular values of the argument (z = 1) was started by 



David Broadhurst [14] in the context of Euler-Zagier sums (or multidimensional zeta val- 
ues) [15]. This activity has received further consideration for another, physically interesting 
point, z = 1/4 (see the relevant Appendix in Ref. ||, |l0|), and also for the "primitive sixth 
roots of unity" (see Ref. [^]). Over time, other types of sums 2 have been analysed in a 
several publications: harmonic sums |, [l7j, generalized harmonic sums [|l^, 12|, binomial 
sums iPH [T^ ] and inverse binomial sums [12, |I9|| . 

The introduction of new functions, such as multiple polylogarithms, (see Appendix [A]) 
independently in mathematics and physics ju], |2^, 22, [2j|, 3 allows us to derive a set of 
universal algorithms for the simplification and construction of the analytical coefficients of 
the Laurent expansion of a large class of hyper geometric functions. (For details, see Refs. 
|H| m, 0, H, ||, |7], ||, ||.) Recently, similar probl ems have also drawn the attention of 



mathematicians [28, 30]. However, the general solution of this problem remains unknown. 

The multiple series representation has further applications in the framework of Feyn- 
man diagram calculations. In particular, the Smirnov-Tausk approach |3| (see also Ref. 
|3l| l) was very productive for constructing the analytical coefficients of the e-expansion (fi- 
nite part mainly) of Feynman diagrams depending on one or two massless (ratio of massive) 
kinematic variables. Presently, there are several computer realizations of this approach ||. 
In the framework of this technique, the Feynman parameter representation |l| for a dia- 
gram is rewritten in terms of multiple Mellin-Barnes (contour integral) representations, 
resulting in expressions for which a Laurent expansion about e = may be constructed 
explicitly, using gamma functions and their derivatives. The results may be summed an- 
alytically or numerically, typically leading to the same sums as in the construction of the 
e-expansion of hyper geometric functions: the (generalized) harmonic sums and (inverse) 
binomial sums. Inverse binomial sums typically arise from massive loops; see Refs. |], |32| . 
Another source of multiple sums in Feynman diagrams comes from the Frobenius series 
solution of a differential equation |33|| . Other classes of sums have been considered as well 4 . 

Analytical results are possible when these sums can be evaluated explicitly. For the 
analysis of (generalized) harmonic sums, the nested sums approach [[l?], [l8| permits the re- 
duction of any type of (generalized) harmonic sum to a set of basis sums. The analytical 
evaluation of these basis sums is an independent problem. (See, for example, Ref. Jl^].) 
The Generating function approach [^] is a universal method for analytically evaluating ar- 



One of the classical tasks in mathematics is to find the full set of parameters and arguments for which 
hypergeometric functions are expressible in terms of algebraic functions. Quantum field theory makes a 
quantum generalisation of this classical task: to find the full set of parameters and arguments so that the all- 
order e-expansion is expressible in terms of known special functions or identify the full set of functions which 
must be invented in order to construct the all-order e-expansion of generalized hypergeometric functions. 
2 See Eq. @ for clarifying of terminology. 

3 Hyperlogarithms have been considered by Kummer, Poincare, and Lappo-Danilevsky; see H|. The 
interrelation between hyperlogarithms and multiple polylogarithms has been discussed in J25[. 

^Finite harmonic sums are another class, on which more details may be found in Ref. fl34]. However, 
there presently does not exist an appropriate generalization of multiple (inverse) binomial sums to finite 
harmonic sums. Some recent attempts in this direction have been discussed in Refs. Jl8|, fj5| . 
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bitrary sums, which was successfully applied (see section (2.3) in Ref. [12]) for an analysis 
of multiple (inverse) binomial sums [|ll], |l^] . The generating function approach allows us to 
convert arbitrary sums to a system of differential equations. The question of the express- 
ibility of the solution to this differential equation in terms of known (special) functions is 
not addressed by this approach. In particular, the partial results of Refs. [0, |ll|, 12, 13] 
were restricted by attempts to express the results of the calculation in terms of only classi- 
cal or Nielsen polylogarithms [37, BMj. It is presently unknown what type of sums (beyond 



generalized harmonic sums) are expressible in terms of known special functions. 5 

The aim of this paper is to prove the following theorems: 
Theorem A 

The multiple (inverse) binomial sums 

oo _. j 

Y,^k Z - s ^-i)---s ap (j-i) 1 (i.i) 

where k = ±1, S a (j) = Y%=i W * s a harmonic series, and c is any integer, are expressible 
in terms of Remiddi- Vermaseren functions with 

(1) for k = 1 : (i) c > 2 rational coefficients; (ii) c < 1 ratios of polynomials; 

(2) for k = — 1 : (i) c > 1 rational coefficients; (ii) c < ratios of polynomials. 
Theorem B 



The all order e-expansion of the generalized hypergeometric functions [41] 

pF^ (A+ae;B+b£, \+h\z} , (1.2a) 
p F p .i (A+ae, l+I 2 ;B + be;z^j , (1.2b) 

where A, B are lists of integers and I\, I 2 are integers, are expressible in terms of the 
harmonic polylogarithms with coefficients that are ratios of polynomials. 

The paper is organised as follows. In section we will prove Theorem A. In section |3|, 
the results of Theorem A will be applied to hypergeometric functions to prove Theorem 
B. Section |] is devoted to a discussion of an algorithm for the reduction and analytical 
evaluation of generalized multiple (inverse) binomial sums. Appendix [A| contains some 
basic information about relevant special functions. 



5 There is not universal agreement on what it means to express a solution in terms of known special 
functions. One reasonable answer has been presented by Kitaev in the Introduction to Ref. where 
he quotes R. Askye's Forward to the book Symmetries and Separation of Variables by W. Miller, Jr.,po[] 
which says "One term which has not been defined so far is 'special function'. My definition is simple, but 
not time invariant. A function is a special function if it occurs often enough so that it gets a name". Kitaev 
adds, "... most of the people who apply them . . . understand, under the notion of special functions, a set of 
functions which can be found in one of the well-known reference books. . . ." To this, we may add "functions 
which can be found in one of the well-known computer algebra systems." 
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2. Analytical evaluation of a basis of multiple (inverse) binomial sums of 
arbitrary weight and depth 

The main purpose of this section is to prove Theorem A. In the first subsection, we will 
consider differential equations satisfied by multiple (inverse) binomial sums, and use the 
analytical properties of such sums to derive two useful lemmas. In the second subsection, 
we prove auxiliary propositions for the separate cases of multiple binomial sums and inverse 
binomial sums, and use them to complete the proof of Theorem A. 

2.1 Some analytical properties of multiple (inverse) binomial sums of arbitrary 
weight and depth 

Let us define the multiple sums 

CO j 

Z ( al.,a p ; l:\ ,— ,b q ;c( " ' = £ ^^(j -!)••• S* p ( j - l)Sfc(2j - 1) • • • S bq (2j - 1) , 

i=i l/J J 

(2.1) 

where S a (j) = Y^k=i & * s a harmonic series and c is any integer. For particular values of 
k, the sums (|2.l|) are called 

( 



k 



generalized harmonic 

1 inverse binomial ) sums 
-1 binomial 



The case S Q1 ... d • q ... o-c( u ) * s called a harmonic sum. The number w — c + a\ + • • • a p + 
b\ + • • • b q is called the weight and d = p + q is called the depth. 

The general properties of multiple sums can be derived from their generating functions. 
Let us rewrite the multiple sum (p. II) in the form S^f- (u) = y^JLi u^rvfl (j) , where a = 



a;b;c 3 a;6;c 

(ai,...,a p ) and 6 = (6i,...,& 9 ) denote the collective lists of indices and rjil (j) is the 
coefficient of u 3 . In order to find the differential equation for generating functions of 
multiple sums it is necessary to find a recurrence relation for the coefficients (j) with 

(k) 

respect to the summation index j. Using the explicit form of rj\l (j), the recurrence 
relation for the coefficients can be written in the form 6 

[2(2j + 1)] \j + ir S<g ;c (j + 1) = f^JJ) + r$0") , (2.2) 



where the "remainder" r\l(j) is given by 

fc p <7 



2 /) = n m - i )+r ar ] x n [^(2i - i)+(2j)- 6 '+(2j+i)- 

J / ' r =l 1=1 

-flf[S ar (j-l)S bl (2j-l). (2.3) 



r=l Z=l 



We would like to point out that Eq. (2.2) is valid for an arbitrary integer fc and c — fc > 0. In the case 
c — fc < 0, the proper term will be generated in the r.h.s. of the equation. 
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Multiplying both sides of Eq. (2^2) by u J , summing from j = 1 to oo, and using the 
fact that any extra power of j corresponds to the derivative u{d/du) leads to the following 
differential equations for the generating functions Eff2 (u) (see Ref. [[0|): 



»4-2] («^) C e£>(«) = 2^,0 + 2 (2u;fc + 1) , (2.4c) 



where R^l(u) = Y^T=i u ^ r<y ^{3) an d ^ s the Kronecker J-function. The boundary con- 



ditions for any of these sums and their derivatives are 



%; ;c (°) = ' j = 0,1,2, ■ 



(2.5) 



From the analysis in Refs. [11, |T^, [T^, 19], we have deduced that the set of equations 
for the generating functions has a simpler form in terms of a new variable. For multiple 
inverse binomial sums, this variable is defined by 

-V« . (!-y) 2 



y 



u 



and for multiple binomial sums, it is defined by 

1 - y/l-Au 



X 



1 + y/1 -Au 



u 



V 



X 



(i + x) 2 



(2.6) 



(2.7) 



Let us consider the differential equation for multiple inverse binomial sums in terms of 
these new variables. The notation S^ fc 2 (y)[(x)] wm De used for a sum defined by Eq. (|2.1j), 



where the variable u is rewritten in terms of variable y[x] defined by Eq. (|2.6|) [(]2.7|)1: 



,(1) 



a:b:c 



( U (y)) = E « c («) 



S^M = («(y)) = E ( "i } (u) 

a:b:c yAJ a:b:c V a;fe;c V ; 



u=u{y) 
u=u(x) 



(21 



In terms of the variable y, equation ( |2.4a| ) may be split into sum of two equations 

(2.9a) 



(2.9b) 



with boundary condition 



E« (1) 



. 



Equation ( p.9a|) could be rewritten as 

l-y d_^ c - j 
l+y V dy 



t& (y) 

a:b:c 



^lXv) 

a;b;j 



(2.10) 
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or in equivalent form 



i ~ y v4-Y j ~ 1 * ( &M= r^f^---)s«.(,). .2.11. 



l + y dy J a;6 ;c J 1 \l — y y J 

From this representation we immediately obtain the following lemma: 
Lemma A (see Ref. Jl9[) 

If for some integer j, the series E^i (u) is expressible in terms of Remiddi-Vermaseren 

S;b;j 

functions ( |A.4b| ) with rational coefficients, then the sums E^g _. + .( u ) for positive integers i 
can also be expressed in terms of functions of this type with rational coefficients. 

In a similar manner, let us rewrite the differential equations for the generating function 
of the multiple binomial sums as 

—x4-Y ^ }) (X) = — ^(x) , (2.12a) 
1-X d x J l-x ^ W ' ^ ' 

\^K^ { ^ s A^A +i ) R ^- (2 - i2b) 

The first equation may be rewritten as 

or in an equivalent form 

^X^r^^Cx) =l X dx(--^A • (2-14) 



1-X S;6;c Jo VX 1 + X/ 



In this case, the boundary condition (2.5) is unchanged, and we can make a statement 
similar to the previous one: 

Lemma B (see Ref. fl 
If for some integer j, the series E^-. (u) is expressible in terms of harmonic poly logarithms 

QA.4b| ) with rational coefficients, then the sums E^,g .(it) for positive integers i can also 
be expressed in terms of harmonic poly logarithms with rational coefficients. 



2.2 Analytical evaluation of multiple (inverse) binomial sums of arbitrary 
weight and depth 



Let us now consider the special case of sums (2.1) including only products of harmonic 
sums, and show that they are expressible in terms of Remiddi-Vermaseren functions 7 with 



7 These sums are related to the multiple sums 



argument (j — 1); see Eq. In agreement with Ref. [12], we will denote such a sum as 

^ii — a p ---m( n )- ^ n this case, the non-homogeneous term r?_(j) of differential equation 
(2.4a) is again expressible in terms of sums of the same type, „ b . _. m (u), but with 
smaller depth: 



2 f) r (k) 



= n - !)+r a i -n - 1) > * = ±i 



(2.15) 



r=l r=l 

We shall start with the case of inverse binomial sums, k = 1: 



CO ^ j 

j=i ^ J 

In order to prove Theorem A for inverse binomial sums, we will prove an auxiliary 
proposition: 

Proposition I 

For c = 1, the inverse binomial sums are expressible in terms of harmonic poly logarithms 
with rational coefficients c rj g times a factor (1 — y)/(l + y): 



^ii,- ,a p ; -;l( u ) 



_ l-y 

u=u(y) 1 + y 



TT^ Yl °r,s^ r y (y) 
1 + y ^ \s) 



(2.16) 



where r + s\ + • • • + s/- = 1 + «i + • • • + dp (weight of l.h.s. = weight of r.h.s.). 

Substituting expression ( p. 16 ) in the r.h.s. of Eq. ( 2,11| ), setting j = 1, and making 
trivial splitting of the denominator, we get the following result: 

Corollary A: 

For c > 2, the inverse binomial sums are expressible in terms of harmonic poly logarithms 
with rational coefficients d r g: 



,(i) 

J a 1 ,--- ,a p ; — ;c 



(tl) 



u=u(y) 



^2 d r ,s^ r V Li(J) (y) , c > 2 



(2.17) 



where r + s\ + • • • + s& = c + a\ + • • • + a p (weight of l.h.s. = weight of r.h.s.). 



Proof: 

Let us consider inverse binomial sums of depth 0: 



. Ju) = ^ 



It was shown in Ref. []|] that for any c > 2 this sum is expressible in terms of generalized log- 
sine functions [37] which could be rewritten [10, ^2| in terms of Nielsen polylogarithms. [38] 
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Here we will present an iterated solution for the case of interest. The system fl2.9|) has the 
form 



1-y d 
l+y V dy 



c-1 



For c = 1, we immediately get the relation 



S^ 1 ? .1 (w) = In y , 



(2.18a) 
(2.18b) 

(2.19) 



which coincides with Proposition I and can be readily transformed into the form of 

1 



Eq. (|2~T0|) : 



1-2/ d V 2 v(!) f ^ 



ln 2 y . 



l + y° dy ; — 2"-- < 2 ' 20 » 

The iterated solution of this differential equation for an arbitrary integer c > 2 is expressible 
in terms of Remiddi-Vermaseren functions with rational coefficients 8 (in accordance with 
Corollary A). 

For sums of depth 1, i.e. 



i (?) ? 



1=1 



the coefficients of the non-homogeneous part are equal to inverse binomial sums of the zero 
depth, ( 2j ) (j) = , and Eqs. (§J) take the form 



1 + V y dv °avAy) 



l+y dyj 



For c = 1 the system of equations ( |2.21|) takes the simplest form 



(2.21a) 
(2.21b) 

(2.22a) 
(2.22b) 



Let us now consider the case a\ = 1. Using Eq. ( 2.19 ), we derive from Eq. ( p. 22 ) 
a[f_ (y) = X - In 2 y - 2 In y ln(l + y) - 2Li 2 (-y) , 

i.e., result is expressible in terms of harmonic polylogarithms. For a\ > 2, the r.h.s. 
of the second equation (2.22b) is expressible in terms of harmonic polylogarithms with 



'Compare with the results of Refs. j| [u| ||. 
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rational coefficients (in accordance with previous considerations), so that crj^_(y) will also 
be expressible in terms of harmonic polylogarithms with rational coefficients. Substituting 
these results in the first equation (2.22a), we obtain results in accordance with Proposition 
I. For c > 2, the desired result follows from Lemma A. 

We may complete the proof by mathematical induction. Let us assume that Propo- 
sition I is valid for multiple inverse binomial sums of depth k: 



ai,— Ak\ 



oo 1 3 

;i («) = E 72^-^0' 

.7=1 Vj) J 



1 + 2/ 



^c r ^ln r y Li/s\ (y) , 



u=u{y) 



(2.23) 



where Li/s\ (z) is a coloured polylogarithm of a square root of unity, s = s%, ■ ■ ■ ,Sk, and 
r + s\ + • • • s p = c + a% + • • • + ajfc. Then for c > 2, Corollary A also holds for multiple 
inverse binomial sums of depth k: 



^Oi,-- ,a fc ;-;c( n ) 



oo ^ 

X, T^YT^ 5,01 G? - " " " ^a* (■? ~~ 



S (?) i 



^c r ^ln r y LL^ (y) , 



u=u{y) 



(2.24) 



For the sum of depth fe+1, the coefficients of the non-homogeneous part may be 
expressed as linear combinations of sums of depth j , j ' = 0, • • • , k, with integer coefficients 
and all possible symmetric distributions of the original indices between terms of the new 
sums: 



i-y ± 

1+y dy 



c-1 



S 



(i) 



d 



k 



l-V (i) 

1 



y 3=1 \ j) P=0(ix,-,iin-l) V 



^ O'-i) ■••^O'-i) 



*p+l+-"»fc+l 



(2.25a) 



, (2.25b) 



where the sum over indices {i\, - - - ik+i) is to be taken over all permutations of the list 
{fli,--- , afe+i). If ip+i + ■••ijfc+i > 2, the r.h.s. of Eq. ( 2.25b| ) is expressible in terms 
of harmonic polylogarithms of weight k with rational coefficients; see Eq. fl2.24|) . As 
the result of integrating this equation, ... 0fc i; _(y) also will be expressible in terms of 
harmonic polylogarithms of weight fc+l with rational coefficients. 

If ip+i + • • • ik+i = 1) the r.h.s. of Eq. ( |2.25b| ) is expressible in terms of harmonic 



polylogarithms of weight k with a common factor (1 — y)/(l+y); see Eq. ( 2.23 ). The result 
of integrating this equation again will be expressible in terms of harmonic polylogarithms 
of weight fc+l with rational coefficients: 



T (1) 
'a\,- 



i a fc+i; 



Xv) 



r^G"ifi)^ c ^ inrtLi (j) (t) - 
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For c = 1, direct substitution of the previous results into ( |2.25a| ) will show that Proposi- 
tion I is valid at weight k-\-l. In this way, the Proposition I is proven for all weights. 
Then for c > 2, Corollary A is also true for multiple inverse binomial sums of depth fc+1. 



... _. c (u), we can derive results for a similar sum with c < l. 9 Thus, Theorem A 
proven for multiple inverse binomial sums 



Applying the differential operator u-^ = — j+^y-^ repeatedly I times to the sum 



is 



10 



Let us now consider the multiple binomial sums 11 , {k = —1), T,^}. ap . _. c (u) 

E( 2 /)^(i-i)-^o--i)- 

In order to prove Theorem A for binomial sums, we will first prove the following 
auxiliary proposition: 

Proposition II 

For c = 0, the binomial sums are expressible in terms of harmonic polylogarithms and have 
the following structure: 



(u) 

oi,- ,Op; -;0V"7 



w=w(x) 



1 



1-X 



In 



r X Lip) (x) 



(2.26) 



where r + si + • • • + s^ = 1 + a\ + • • • + a p (weight of l.h.s. = weight of r.h.s.) and c r 

and d r s are rational numbers. 



Substituting the expression fl2.26|) in the r.h.s. of Eq. (2.14) and setting j = 0, we get 
Corollary B 

For c > 1, the binomial sums are expressible in terms of harmonic polylogarithms with 
rational coefficients d r g: 



u=u(x) 



^2 d r ,gln r x LiQ (x) , c > 1 



(2.27) 



where r + s\ + • • • + Sk = c + a\ + • • • + a p (weight of l.h.s. is equal to weight of r.h.s.). 



9 Some particular cases of sums of this type were considered also in Ref. 

10 All multiple inverse binomial sums up to weight 4 were calculated in ref. see Table I in Appendix 

C. 

11 These sums are related to the multiple sums 



PI- 

ef~jl|]; 



/ '111 | 



ni >7i2 >■ '-np — 1 
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We start again from the multiple binomial sums of depth 0, 



2 A u 



J J J 



In this case, Eqs. ( |2.12 ) have the form 



where the factor may be written as 

1 + X 



For c = 0, we obtain 



i-x i-x 



i-x 



1 



(2.28a) 



(2.29) 



(2.30) 



which coincides with Proposition II. Substituting this result into r.h.s. of Eq. ( 2.14| ) we 
find 



(2.31) 



The results of iterated integration, for c > 1 and boundary condition defined by Eq. fl2.5|), 
are expressible in terms of generalized poly logarithms (|A.2|) with rational coefficients (see 
Corollary B). For the sums of depth 1, 



we have 



1 + X ±_ 

1-X dx 



,(-1) 



i=i 



7/7 — * 

J 7 ^ i=l 



^(l+x) 2 ^^(x)=2S^ ) ;ai _ 1 (x)+S^ ) ;ai (x) 
Integrating by part in Eq. (2.32b) and taking into account that 12 



(2.32a) 
(2.32b) 



dx 



x i + x 



12 This relation follows from the differential relation 

d_ 
du 



U ^. E W („) = E M («) 
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we obtain 



x dt„(_i) 



(2.33) 



Using this results in the r.h.s. of Eq. ( 2.32a| ), we have 

Let us set c = 0. It is necessary to consider two cases: (i) a\ = 1 and (ii) a\ > 2. For 



a\ = 1, we can use the explicit results ( [2.30D and ( 2.31 ) to get 

SQW = 21n(l - X ) " 21n(l + X ) ~ ~7~~l_ Ml " x) , 

in accordance with Proposition II. For ai > 2 the r.h.s. of Eq. ( |2.32b| ) is expressible 
in terms of harmonic polylogarithms with rational coefficients, so that Eq. (2.34) is also 
expressible in terms of harmonic polylogarithms with rational coefficients in accordance 
with Proposition II. 

For c > 1, the desired result follows from Lemma B: 



1+X d 



c-l 



In particular, for a% = 1 we have 

^X^J 1 <_ 1 ; ) c (x)=2Li 2 (-x) + 21n 2 (l + X ) + 2Li 2 (x) . 



o h Jo h 



i(* 2 ) • (2-35) 



(2.36) 



Let us assume Proposition II is valid for multiple binomial sums of depth k, and 
prove the proposition for depth fc+1. Thus, we assume that 



) = E(^V s ».w- 1 )--- s »»G-i) 



u=u(x) 



i-x 



Cp 5 s "I - C^p,s 



ln P xLi(<?) (x) , 



(2.37) 



where (x) is a coloured polylogarithm of a square root of unity, s = (si, • • • , Sfc), and 
p + si + - ■ ■ s p = a%-\ — - + afc. Then for c > 1, Corollary B also holds for multiple binomial 
sums of depth k: 



2 A « J 



3/3 



S ai (j ~ I) ■ ■ ■ S a Jj ~ 1 



u=u(x) 



(2.38) 



For a sum of depth fe+1, the coefficients of the non-homogeneous part are expressed 
as linear combinations of sums of depth j, j = 0, • ■ • , k, with an integer coefficients and 
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all possible distributions of the original indices between terms of new sums, multiplied by 
a factor (2j + 1): 



1+X d 



x— Vs M) ( x) = l±^ CJ (- 1 ) ( X ) 



i-x rf x. 

^( 1 + X) 2 ^^ 1, ,a fc+1 ;-(X)=E^ + 1 : 



(2.39a) 



E E 

P=0 (i 1( -,i k+ i) 



5 ix (j-1)-- -S ip (j-l) 



(2.39b) 



where the sum over indices (i\, ■ ■ ■ ,ik+i) is to be taken over all permutations of the list 

Let us denote the sub-list of length p as I = (ii, • • • , i p ) and define the sum of the 
remaining indices as J = i p+ i + • • • + ik+i, so that the second equation ( 2.39b ) can be 
written as 

^(i+x) 2 ^^, afc+i; -(x) = E K^W + CUW 

Integrating by parts, we find 

(x) = T 

I, J 



LJ 



l-X v (-i) 

i + x' 



Substituting this result into the r.h.s. of Eq. (2.39a) we have 

1+X v d VyH) y (-i) M , j+X f Xdt i^(-i) t± \ 



I, J 



(2.40) 



Let us set c = and consider two cases: (i) J = 1 and (ii) J > 2. For J = 1, the 
first term of the r.h.s. of Eq. ( 2.40 ) is expressible in terms of harmonic polylogarithms with 
rational coefficients. The last term of the r.h.s. of Eq. ( [2.40 ) has the structure of Eq. ( 2.37| ) 
so that after integration, it will again be expressible in terms of harmonic polylogarithms 
of weight fc+1. For J > 2, both terms of the r.h.s. of Eq. ( |2.40| ) are expressible in terms 
of harmonic polylogarithms of weight fc+1 (see Eq. ( 2.40| )). In this way, the Proposition 
II is found to be valid at the weight fc+1. Consequently, Proposition II is proven for 
all weights. Therefore, for c > 1, Corollary B is also valid for the multiple binomial sums 
of weight fc+1. 

Applying the differential operator u-^ = jz^x~^r repeatedly I times to the sum 



~(x)> we can derive results for similar sums with c < 0. Thus, Theorem A 



is proven for multiple binomial sums. 



13 



13 All multiple binomial sums up to weight 3 were calculated in ref. jll], [ill; see the proper Appendixes. 



- 13 - 



3. All-order e-expansion of hypergeometric functions with one half-integer 
value of the parameters via multiple (inverse) binomial sums 

In this section, we turn our attention to the proof of Theorem B. It is well known that 
any function p F p -i(a + m; 6 + k; z) is expressible in terms of p other functions of the same 
type: 

v 

R p+1 (a,b, z) p F p ^i(d + rh;b + k; z) = R k (a,b, z) p F p ^ 1 (a + el; 6 + E k ; z) , (3.1) 

k=i 

where rh, k, e k , and E k are lists of integers and R k are polynomials in parameters a,b, 
and z. Systematic methods for solving this problem were elaborated in Refs. [44, 45]. 
For generalized hypergeometric functions of Theorem B, let us choose as basis functions 
arbitrary p- functions from the following set: 

• for Eq. (1.2a) there are p 2 functions of the proper type: 



pFp—i 



l{l + aie}P- L -\ {2 + diE} L 



{l + e i e}P-Q-\{2 + c i £}^ 
for Eq. ( 1.2b|) there are p 2 — 1 functions of the proper type: 



pFp-i 



{1 + a i£ }P- L , {2 + d,e} L 
\,{l + e i ey-Q- 2 ,{2 + c i eY 



In the framework of the approach developed in Refs. || |l0|, 11, [T^, [T^] , the study of the 
e-expansion of basis hypergeometric functions has been reduced to the study of multiple 
(inverse) binomial sums. It is easy to get the following representations: 



pFp—x 



pFp-i 



{l + ai e} K , {2 + d l e} 1 
\,{l + e i e} R ,{2 + c l e} 



i n Q =1 (i + c sg ) ^ i 
2z nf =1 (i + d i£ )^(^)i 



(4z) 



2j\ jK-R-1 



A 



§,{l+Oi£} A , {2+die} 
{l+eie} R ,{2+ Cl e}Q 



2 lff =1 (l+c. e ) ^ 
zD? = i(l+dis)^ 



V 



3)3 



;K-R-1 ' 



(3.2a) 
(3.2b) 



where the superscripts K, L, R, Q show the lengths of the parameter lists, 



A = exp 



E 

k=l 



w k j +S k (n-l)t k 



\ [ j- +t 1 S 1 (n-l) ) +OC-) . (:U) 



S a (n) = X]j=i 1/ ' 3 a ls a harmonic sum, and the constants are defined as 

tk = Cfc + Ek — Ak — D k , Wk = C k — D k , 

where the summations extend over all possible values of the parameters in Eqs. ( |3.2|) . In 
this way, the e-expansions of the basis functions (|3.2| ) are expressible in terms of multiple 
(inverse) binomial sums studied in Sect. ||. But all these are are expressible in terms of 
harmonic polylogarithms. Thus, Theorem B is proven. 
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4. Generalized multiple (inverse) binomial sums via derivatives of 
generalized hypergeometric functions 



In physical applications, in particular, within Smirnov-Tausk approach, more general sums, 
in addition to the ones defined in Eq. (1.1), may be generated: 



E 

3=1 



(j + ci)!(j + c 2 ) 
(2j + c 3 )! 



(nj + c 4 ) 



: S ai (mij + bi)---S ak (m k j + b k ) 



where {a*}, {bj}, {c k }, {m^}, n are integers and k = ±1. The procedure of finding the 
proper differential equation (see Refs. [12, 36] for a detailed discussion) can be applied 



to analytically evaluate any of these new sums. Another approach is based on extension 
of the algorithm of nested sums 1 17, |l8| for the study of the algebraic relations between 
these sums. However, there is a third approach arising from the possibility of reducing an 
arbitrary generalized hypergeometric function to a set of basis functions with the help of 
the Zeilberger-Takayama algorithm described by Eq. (3.1). 

To be more specific, let us divide both sides of Eq. ([O]) by i? p +i(aj, bj, z) and construct 
the e-expansion for the hypergeometric functions described in Theorem B. The r.h.s. of 
this relation is expressible analytically in terms of harmonic polylogarithms with polynomial 
coefficients. The l.h.s. can be used as a generating function for generalized multiple (inverse) 
binomial sums. Using a standard form for the Taylor expansion of the Gamma function, 1 



[m 



ae) 



exp 



E 

k=l 



-ae) 



[S k (m+j-l)-S k (m-l)] 



where fa) 



p+xFp 



r(a + j)/T(a) is the Pochhammer symbol, we obtain 

nf =1 (m, + a l£ )j ^=i~ L {Pi + h)j 



{m l+ a l e} L ,{p i + ±} p + 1 - L 



{n k +b k e} K ,{ qi + ±} p - K 



E 

3=0 



E 



II/' :!"'/*, Vr' (2/', + n 



P-K 



where the mi,n k ,pi,qj are integers and 



J (2l s +l) 2j 



A 



(4.1) 



A 



exp 



K 



E 

.k=l 



[a^Skimi+j-^-afShirrii-l) 



i=l 



14 The relation between harmonic sums S a (j) and derivatives of the function ip(z) = 4-\nT(z) is 

= (-l) fc (fc - 1)! [C* - S„(j - 1)] , fc > 1. 
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Setting K = L = P in Eq. (|4.l|), we get generating functions for generalized multiple 
binomial sums: the derivatives 



n 

Lk 



d_ 

dai 



d_ 

dh 



p+i 



F P 



{bkY 



(4.2a) 



ai=mi;b k =n k 



lead to terms in the epsilon expansion of the form 



(4.2b) 



where the Im are integers from the lists {mi} L and {nk} K . For L = P + 1 and K = P — 1 
we get generating functions for generalized multiple inverse binomial sums: 



n 

Lk 



dai 



_d_ 

dbl 



p+iFp 



{ai} P+1 
{b k } p -\q+\ 



ai =rai $^=71.^ 



~ (9 + 1), (Azy nl+Hm^ tt 



(4.3) 



M=l 



For K = P and L = P + 1 we get generating functions for generalized multiple harmonic 
sums: 



nt-YY- 



db k 



p+iFp 



({ai} P+1 




{ {h} P 


■) 



ai=mi;b k =n k 



00 i TtP+Xt \ 

^o- 7 " u *=iWj m=1 



(4.4) 



Instead of one hyper geometric function, we could consider a linear combination of 
the functions of the same type. Such a combination is also reducible and expressible in 
terms of our basis functions. Combining the proper set of hypergeometric functions, we 
could expect that any individual sums, 15 of the type described by r.h.s. of Eqs. ( [I~2j ) - 
Q4.4| ) are expressible in terms of generalized (harmonic) polylogarithms with polynomial 
coefficients. 16 



15 Using the results of the all-order e-expansion for Gauss hypergeometric functions [ [l3[ pj| we could 
consider a series of type (2.1). 

16 In particular, all sums presented in Ref. [[|6| are reducible in terms of our basis sums or sums studied 
in Ref. ]l2| ]. Indeed, taking into account that 

n + 1 / 2n + 2 s 



(2n + l)( 2n 

V n 

and shifting the index of summation we have 

CO n oo 



2 V n + 1 



' ^X 3 (j-l)F ? (2j-l) - X 3 (0)Y S (1) 



(?) 



(4.5) 
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These arguments suggest a criterion for what type of generalized multiple (inverse) bi- 
nomial sum are expressible in terms of harmonic polylogarithms with coefficients that are 
ratios of polynomials. This is just the beginning of a general analysis, but the corresponding 
analysis for harmonic sums is already known to be valid. |l7j Unfortunately, existing com- 
puter algebra algorithms [47| do not allow us to identify the multiple series with derivatives 
of hyper geometric functions or their combinations. It is still matter of personal experience, 
but this approach looks very promising and is worthy of further analysis. 



5. Discussion and Conclusions 



We have constructed an iterative solution for multiple (inverse) binomial sums defined by 
Eq. ((Ll]). It was shown that by the appropriate change of variables, defined by Eqs. ( |2.6| ) 
and ( |2.7D , the multiple (inverse) binomial sums are converted into harmonic polylogarithms 
(see Theorem A). Symbolically, this may be expressed as 



1 u? 

3=1 \j) J 
00 1 3 
3=1 \j > J 



1 — y 

rqr - $Z c p,^ lnP y L \i) (y) > ( 5 - la ) 



u=u{y) 



p,s 



= ^SpjhPy (y) , c > 2 (5.1b) 

u=u(y) p,s 



and 



E 

i=i 

oo 

E 

3=1 



2j 
3 

2j\ v? 



3/3 



uPS ai ti-l)---S a „V-l) 
S ai (j-l)...S ak (j-l) 



u=u( X ) P> g 



+ d 



p.s 



hl p 



X Li(j) (x) , (5.2a) 



= J2 ~ c p,s^ p X Li(J) (X) , c > 1 (5.2b) 

u=u(x) P ' S 



where c is a positive integer, c p g, c p g and d p g are rational coefficients, the weight of 
l.h.s. = weight of r.h.s., Li^ (x) is the coloured multiple polylogarithm of a square root 



where Xs(n) = H 1 k - 1 S ak {n) and Vs(2n + 1) = IT k=1 Sa k (2n + 1), are products of harmonic sums, with the 
vector a having r components. As a consequence, Xg(0) = and Yj(l) = 1. In this way, any sums described 
by Eq. (4.5) may be reduced to sums of type (2.1), and for YAj) = 1 they are reduced to the sums studied 
in the present paper. Another possible generalization of the sums considered here is 



Due to the depth reduction relation, 

r 

P=0 (ii,— ,ik+i) 



2 ™ ] ~> 

X 3 (n + l)Y g (2n+l) = _^ 7 -^_X s C7)n( 2 i-l)-- x '3(lW) ■ ( 4 - 6 ) 



3=1 



1 S n (j-i)...g ip (j-i) 
p!(r — p)! j'p+i+'-'t+i 



sums of type (Hji) are also expressible in terms sums of type (4.5). 
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of unity, 



I" 

i=l 



is a harmonic series. The mappings (5.1), (|5.2| ) are defined in the radius of convergence of 
the l.h.s.: 

{4, inverse binomial 
| , binomial 

Unfortunately, one of the unsolved problem is the completeness of the representation 
(|5.l| ), fl5.2|) . In other words, is it possible to express all harmonic polylogarithms in terms 
of multiple (inverse) binomial sums? If not, what kind of sums must be added to get a 
complete basis? Another problem beyond our present considerations is to find the algebraic 
relations among the sums. 

From representation (|5.1|), (|5.2|) , it is evident that some (or all, if the basis is com- 
plete) of the alternating or non-alternating 17 multiple Euler-Zagier sums (or multiple zeta 
values) |15[] ) can be written in terms of multiple (inverse) binomial sums of special values of 
arguments. Two arguments where such a representation is possible are trivially obtained 
by setting the arguments of the harmonic polylogarithms y,x to dbl: 

u = 4, y = -l, (5.5) 

u = - , x = 1 • (5-6) 

Another such point 18 

i 3- V5 

u = -l, y = — - — (5.7) 

has been discussed intensively in the context of Apery-like expressions for Riemann zeta 
functions (see [48] and References therein). For two other points 

u = 1 , y = exp (i^j , (5.8) 
u = 2, y = i, (5.9) 

the relation between multiple inverse binomial sums and multiple zeta values was analysed 
mainly by the method of experimental mathematics. Some of the relations are presented 
in Ref. J5(| and in the appendix of Ref. |jlO| . 

Let us make a few comments about harmonic polylogarithms of a complex argument. 
For the case < u < 4, the variable y defined in (|2.6|) belongs to a complex unit circle, 



r Let us recall that multiple Euler-Zagier sums are defined as 



k 



71 ■ 

ni>n 2 >...>n t >0 j=l 3 



where <Tj = ±1 and Sj > 0. a — 1 is called non-alternating and a — — 1 is alternating sums, correspondingly. 

18 We are thankful to Andrei Davydychev for information about the relation between this point and the 
"golden ratio" , Q , 



3-V5 _ ( l-VS } 2 



2 
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y = exp(iO). In this case, the coloured polylogarithms of a square root of unity can be split 
into real and imaginary parts as in the case of classical polylogarithms. |37| At the present, 
there is no commonly accepted notation for the new functions generated by such splitting. 
In Ref. [50], the multiple Glaishers and multiple Clausen functions were introduced as 
the real and imaginary parts of generalized polylogarithms of complex unit argument. In 
Ref. [10, 27, the splitting of Nielsen polylogarithms was analysed in detail. In this 
case, the real and imaginary parts are reduced to classical Clausen functions, Clj (9) and 
generalized log-sine functions Ls^- (8). Ref. [51] attempts to classify new functions on the 
basis of new LsLsc^j^ (#)-functions. 

In Appendix A of Ref. [12], the iterated representation for Remiddi-Vermaseren func- 
tions of complex unit was constructed. It was observed |8|, |l0|, 12, 



that the physically 

interesting case, representing single-scale diagrams with with two massive particle cuts, 
corresponds to Remiddi-Vermaseren functions ( A. 9 ) with argument equal to a primitive 
"sixth root of unity" , y = exp (zf ) • This gives an explanation of the proper "basis of tran- 
scendental constants" constructed in Refs. [j52j | and [10|, and its difference from the proper 
basis of Broadhurst p6|. Of course, for numerical evaluation of harmonic polylogarithms 



of complex argument, only a series representation is necessary. [p3|| 

Using the results of Theorem A, we have proved Theorems B about the all-order 
e-expansion of a special class of hypergeometric functions. The proof includes two steps: 
(i) the algebraic reduction of generalized hypergeometric functions of the type specified 
in Theorems B to basic functions and (ii) the algorithms for calculating the analytical 
coefficients of the e-expansion of basic hypergeometric functions. The implementation 
of step (i) - the reduction algorithm - is based on general considerations performed in 
Refs. p4| , |4q ]. In step (ii), the algorithm is based on series representation of the basis 
hypergeometric functions defined by Eq. Q4.1| ), The coefficients of the e-expansion are 
expressible in terms of multiple (inverse) binomial sums analyzed in Theorem A. 

Exploring the opportunity to reduce an arbitrary generalized hypergeometric function 
to a set of basis functions with the help of the Zeilberger-Takayama algorithm, we have 
presented in section |] some arguments about one possible generalization of (inverse) bino- 
mial sums (see Eq. (4.4)) which would be expressible in terms of harmonic polylogarithms 
with coefficients that are ratios of polynomials. 
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Appendix: 

A. Zoo of special functions 



For completeness, we will present the definition of a set of new functions, such as 
multiple poly logarithms 19 



z m-i z m 2 _ . , 

U klM> ... ikn (z 1 ,z 2 ,--- ,z n )= Yl "Vfe \- ( A -!) 

rn 1 >m 2 >-rn n >0 m l m 2 ' " m ™ 

Special cases of multiple polylogarithms 20 include generalized polylogarithms, defined by 



U klMt ... tkn (z) = }2 -1^-2 ST' ( A ' 2 ) 

. ~^ ^ n rrii rrin ■■■m n 

and coloured polylogarithms of a square root of unity, defined by 21 

. ^-"11 . . . n m n 

Li (?) (z) ^ Li (CT1 , CT2 ,.. , CTfc) (*) = £ 1 " , (A.3) 



mi >m2 > • • -m n >0 

where s = (si, ■ ■ ■ s n ) and a = (<r 1; • • • , <j n ) are multi-indices and cr^ is a square root of 
unity, a k = ±1. The extension of coloured polylogarithms of square root of unity ( |A.3| ) 
by inclusion of powers of logarithms, ln fc z, leads to harmonic polylogarithms or Remiddi- 
Vermaseren polylogarithms (or functions) |22fl . These can written in the following form: 



h a( z ) = ]>2 c P ,k lnPz Li *i,fc2,-,fc™ (*) > ( A - 4a ) 



= ^ c^ln* z Li (?) (z) . (A.4b) 



where in the first equation, ( |A.4a ), the vector A includes only components and 1, and in 
the second, ( A.4b| ), —1 components are included. The coefficients and eg are rational 



numbers. In Eq. ( A. 4 ) the weight of the l.h.s. = the weight of the r.h.s. 

Recall that generalized polylogarithms ( |A.2[) can be expressed as iterated integrals of 
the form 

, . f z dt dt dt dt dt dt dt dt . . . 

Lit.. ... u (z) = / — o — o • • • o — o o • • • o — o — o • • • o — o , (A. 5) 

u ' n K ' 7 n t t t 1 - 1 t t t 1 - 1 y ' 

fcx— l times k n — l times 

where, by definition 

dt dt dt dt _ f z dti [ ll dt 2 f tk - 2 dt kl -i f^ 1 ' 1 dt kl 



o — o • • • o — o- 

o t t 1 — t Jq t\ J t% Jq t kl -i Jq 1 — t kl 



(A.6) 



fei-i times 



19 For a review, we recommended Ref. |54[ . 

20 Our notations corresponds to Waldschmitd's paper of Ref. IB3 ]. 

21 We call n depth, and k = ki + fa + • • • + k n (s = si + S2 + • • • + s n ) the weight. 
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The integral ( |A,5| ) is an iterated Chen integral [55| w.r.t. the differential forms ujq = dz/z 



and oj\ = tt;, so that 




J o 

The coloured polylogarithms (Eq. |A.3| )) also have an iterated integral representation w.r.t. 
three differential forms, 



dy 

1 

y 

ady 



a = 



a = ±l 



(A.8) 



1 - ay 



so that 




i 



Li 



f CT l.°"2 



•-.'k) (y) 



,*2-l 





(A.9) 
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